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We investigate the superconductivity of a three-dimensional d-p model with a multilayer 
perovskite structure on the basis of the second-order perturabation theory within the weak 
coupling framework. Our model has been designed with multilayer high-Tc superconducting 
cuprates in mind. In our model, multiple Fermi surfaces appear, and the component of a 
superconducting gap function develops on each band. We have found that the multilayer 
structure can stabilize the superconductivity in a wide doping range. 
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1. Introduction 

In the last decade, high-Tc superconducting cuprates (HTSCs) with multilayer structure 
have been investigated extensively by various experimental methods. ^"'^ Multilayer cuprates 
typically have higher Tc than single-layer ones. The nuclear magnetic resonance (NMR) mea- 
surement systematically reveals the local hole concentration on each layer of the multilayer 
HTSCs. In previous excellent studies,^' ^ each layer has been crystallographically classified into 
an outer or inner Cu02 plane in a unit cell. The authors found the relationship between the 
difference in the local hole concentration among these two types of planes and Tc, and they 
considered the condition in which Tc can be maximized. The results suggest that multiple 
Fermi surfaces are involved in superconductivity, and that an extensive theory on multilayer 
materials should be constructed on the basis of the model with multibands. 

The above NMR study mainly revealed the characteristic of the multilayer compounds for 
n > 3, where n is the number of Cu02 planes per unit cell. The characteristic feature of multi- 
band systems can also appear explicitly in a bilayer compound. In Bi2Sr2CaCu208-i-<5 (Bi2212), 
which is another typical bilayer material, the high-resolution angle-resolved-photoemission 
spectroscopy (ARPES) successfully revealed the doubling of a band near the Fermi level. ^'^ 
This splitting is negligible along the (0, 0) — > (vr, vr) nodal line and maximum at (vr, 0) in 
momentum space. This momentum dependence of energy splitting is qualitatively consistent 
with the LDA prediction for YBa2Cu307,^ which is another bilayer material. This LDA cal- 
culation predicted that a Cu 4s orbital has a transfer integral between that in the other layer. 
This interlayer coupling causes a single band to split into antibonding and bonding bands. 
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A multilayer model has been theoretically studied since the early 1990s by fluctuation 
exchange (FLEX) approximation^ and by quantum Monte Carlo (QMC) simulation^' on 
the basis of the two-dimensional (2D) multilayer Hubbard model. These studies have shown 
that the (i^2_j^2-wave pairing correlations are reduced by inter layer transfer, which is inde- 
pendent of the in-plane momenta, kx and ky. However, when we introduce inter layer transfer 
into our model, it is important to consider the symmetry of interlayer transfer integrals in 
which Cu 4s electrons participate. Liechtenstein et al. considered this point and studied the 
extended Hubbard model with anisotropic interlayer hopping, using the FLEX approxima- 
tion.^^ Although they could not reproduce the experimental result of bilayer materials having 
higher than single-layer ones, their work should be appreciated as the first theoretical anal- 
ysis of the superconductivity of realistic multilayer systems. All the theoretical works have 
been done on the basis of the 2D model Hamiltonian. We feel that we should investigate the 
three-dimensional model Hamiltonian with anisotropic interlayer hopping to estimate of 
multilayer materials. 

In this study, we investigate the superconductivity of multilayer perovskite. We adopt the 
three-dimensional (3D) d-p model with anisotropic interlayer transfers as our model Hamil- 
tonian. In our model, we introduce such a small on-site Coulomb interaction that the second- 
order perturbation theory (SOFT) can be justified. We can treat the superconductivity within 
the weak coupling analysis because, in our model, the effective interaction for Cooper pairing 
is so small that only the electrons on the Fermi surface are involved in the superconductivity. 
The weak coupling formalism for the repulsive interaction model, since the pioneering work 
by Kohn and Luttingcr,^'^ has been developed by many theorists }'''''^'^ Recently, this formal- 
ism was applied, by Kondo, to the 2D Hubbard model with the formulation applicable even 
for the case with a very small effective interaction.-*^^ We apply Kondo's formulation to our 
model with multiple Fermi surfaces, and clarify how the superconducting gap depends on n 
of layers. We can show that the calculation on the basis of 3D model Hamiltonian is requisite 
for the true estimation of of multilayer materials. Our obtained results are not only to 
be compared with the actual of multilayer materials but to be considered as a guide for 
designing materials with high Tc- 

2. Formulation 

We can decompose our 3D d-p model with the n-layer perovskite structure into several 
parts as follows: 

n 

1=1 

- ^JZ^dliJkla+PklaPkla+PklaPkla) 1 ' 
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(1) 

where dki^ (d^;^), p%i„ iPkla) Pkla (Pkla) annihilation (creation) operators for 

d-, p^- and p^-electrons of momentum k and spin a = {t,i} on the Z-th layer, respectively. 
We define n + 1 = 1, and the chemical potential is represented by /j.. The noninteracting parts 
in eq. (1), i.e., and -fff+i i , are represented by 
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In eqs. (2) and (3) we take the lattice constant of the square lattice formed of Cu sites, 
a, as the unit length. Using it, we can represent = 4tpp sin ^ sin ^ , = 2itdpSin^, 
(I = 2itdpsm^, d = t^cos^cos^e'''^'''', and Ct = ^(cos/c^ - cos kyf e^'^''^^ , where = 
c +(ra-i)ci, '''' ~ c +(n-i)ct, ' '^'^ ^^'^ '^^ represent the distances among Cu02 planes, as shown 
in Fig. 1. Hereafter, we use 'IP' and 'OP' as abbreviations for inner Cu02 planes {\ <l < n) 
and outer Cu02 planes (Z = 1 or n), respectively. A; is the hybridization gap energy on the 
Z-th layer between d- and p-orbitals. When n > 3, our model has several crystallographically 
inhomogeneous Cu02 planes, and A; varies according to I. We can define A; as 

A, = | ^ = '°^" , (4) 

where a is the constant chosen so as to reproduce the difference in doped carriers between 
OP and IP. Wc consider only the on-site Coulomb repulsion among d-electrons. Thus, the 
interacting part H[ in eq. (1) is described as 



fefe'q 



(5) 



In eq. (5), N is the number of fc-space lattice points in the first Brillouin zone (FBZ), which 
is equal to the number of Cu sites in the real space. 

In the following part, we assume that only the electrons on the Fermi surface of the same 
band can have singlet pair instability. For our n-layer model, n d-like bands always intersect 
with the Fermi level. Thus, according to the BCS theory, we have the following self-consistent 
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Fig. 1. Crystal structure modeled after trilayered perovskite. Filled (Open) circles represent copper 
(oxygen) cites. Dotted (Dot-dashed) arrows conecting copper cites correspond to (C|) in the 
text. 

equation for the pair function on the A-th d-like band, ^kx- 

^'^A = E ^^.e^ + fe')^^M%EL= (6) 

where i,j = 1, . . . ,n (layer indices) and X,v = 1, . . . ,n (d-like band indices). Vij{q) represents 
the effective singlet pair scattering between a d-electron on the i-th layer and one on the j'-th 
layer, e^u represents the energy dispersion of the z^-th d-like band, and Zi\{k) represents the 
matrix element of unitary transformation. They are obtained by solving the eigenequation for 
the noninteracting part in eq. (1). We set = Agc-^'fcA, where Age denotes the magnitude 
of and represents its fc-dependence on the A-th d-like band. On the basis of Kondo's 
argument, retaining only the divergent term, we can rewrite eq. (6) as 

*feA = 



ijk'v 



(7) 



for very small Age- Equation (7) is a homogeneous integral equation for ^j^x with the eigenvalue 
of 1/ log Age. We are interested in obtaining the most stable superconducting state, thus we 
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must find the eigenvector with the smallest eigenvalue 1/ log Age using eq. (7) when Age 
is maximum. In our previous paper, we confirmed that the most stable pairing state near 
half-filling is the d^2_y2-wave, by a similar approach based on 2D d-p model.^^ Hence, when 
we assume that 

*feA = aA(fc^)(cos kx - cos ky), (8) 

we can safely reduce our original eigenvalue problem for '^kX to an eigenvalue problem for 
ax{kz) in order to seek only the most stable pairing state. Furthermore, considering the sym- 
metry of in eq. (8), we can take 

Vij{q) = U\^j{q) 

t2 f-t _ f . a f. (9) 



within SOFT. In eq. (9), 



fkr] — o 



l-tanh(^^ 
2T 



(10) 



and T denotes the temperature. 



3. Results and Discussion 

In our present analyses, all e^i, and Zix{k) in eq. (7) are first calculated for 64^ A;— points 
on an equally spaced mesh in FBZ for each band. Then, we calculate Vij{k + k') in eq. (7) 
only for k— and fe'— points satisfying the conditions e^x = A* ^k'u = respectively. When 
we calculate Vij{k + k') according to eqs. (9) and (10), we set the temperature T = 0.001 eV~ 
10 K, at which our system can be considered to behave similarly to that in the ground state. 
These calculations have been performed at U = 0.5 eV, where magnetic instabilities cannot 
occur. We take Ca = 0.3 and = 0.7 for all n. Other common parameters are summarized in 
Table I. In order to solve eq. (7) practically, we substitute into both sides of eq. (7) using 



^dp tpp tdd i± ^dp 

1.00 eV -300 meV lOmeV 5meV 2.10 eV 
Table I. Transfer and hybridization gap energies. 



eq. (8) and integrate for kx-, k'^., ky, and ky. Thus, we reduce eq. (7) to the eigenequation for 
ax{kz)- When we solve it numerically by the standard method, we can finally obtain both the 
eigenvalue, 1/logg Age, and the eigenvector, ax{kz)- 

First, we summarize our results on Age vs (5h((5e) in Fig. 2. As discussed in our previous 
paper on the two-dimensional (2D) d-p model, the existence of a Van Hove singularity (VHS) 
at the Fermi level causes a high density of states (DOS) and enhances Age- In the 2D d-p 
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Fig. 2. (a) log Age vs (5h (hole-doped) or (electron-doped) in the cases oi n = 1,2,3,4, and 5. (b) 
Magnification of (a) in the region of 0.0 < (5h < 0.2. 



model the parts of the Fermi surface with VHS are distributed as lines. Therefore Age varies 
drastically in the neighborhood of the doping point at which the Fermi surface has VHS. This 
is in contrast with the case in the 3D d-p model. Although the energy dispersion along the 
c-axis introduced in our analyses is very weak, as indicated in Table I, the parts of the Fermi 
surface with VHS are distributed as points. Thus, the transition of Age in the 3D model is 
milder than that in the 2D model, as seen in Fig. 2. 

Furthermore, in Fig. 2 we can clearly recognize that the enhanced Age prevails in a wider 
doping region with larger n. This result is caused by the multilayering effect introduced as 
described by eq. (4). In order to explain how the multilayering effect occurs, we show the 
eigensolutions, ax{kz), defined by eq.(8), and the Fermi surfaces for the cases with n = 3,4, 
and 5 in Figs. 3, 4, and 5, respectively. 

In Fig. 3 the highest amplitude of the eigensolution appears in the band with A = 2, 
whose Fermi surface most closely approaches the VHS points, i.e., (±7r, 0, [— vr, vr]) and 
(0, ±7r, [— TT, vr]). Thus, the band with A = 2 has the largest DOS near VHS points and is 
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Fig. 3. The results for n = 3 and = 0.127. (a) The eigensolution a\{kz). (b) The Fermi surface 
projected onto the plane with kz = 0. 



dominant in the superconductivity. For the same reason, the largest amplitude of the eigen- 
solution appears in the band with A = 3 and in the one with A = 4, as shown in Figs. 4 and 
5, respectively. 

When we change the amount of doped carriers, the Fermi surface should be transformed. 
As a result, another band could then have the largest DOS near the VHS points and dominate 
the superconductivity. This possibility should be further increased more if our model has more 
alternative bands. Hence, the enhanced Asc tends to prevail in a wider doping region with 
larger n. This tendency should remain when n becomes much larger, as long as the conventional 
Fermi surface can be defined. However, the largest value of Asc should be saturated toward 
the intrinsic value for n ^ oo. 

Hereafter, we turn our attention to the maximum of the n-layered materials, which 
would be proportional to the maximum Age in our calculated results. In several real materials, 
the largest is achieved when n = 3 or n = 4, and is rather low when n = 5.^'^ For such 
materials our assumption that the well-defined Fermi surface exists might not be valid. For 
example, in the five-layered compound HgBa2Ca4Cu50y, the inner Cu02 planes turn out to 
be antiferromagnetic on account of the strong electronic correlation.^'^ Concerning the strong 
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electronic correlation, other theoretical works based on the 2D multilayer t-J model have been 
extensively carried out by Mori et al.^^''^^ Their approach would be better for explaining the 
results for such compounds. 

Although our results on Tc are not consistent with those for several real materials, our 
conclusion on the multilayering effect is clearly applicable to other real materials. Indeed, 
(Cu,C)Ba2Ca3Cu40i2+j/ (Cul234), has a high Tc even though it is in the heavily overdoped 
region. Cul234 has been revealed, by NMR experiment, to have doped holes that are almost 
uniformly distributed into each layer. ^'^ Thus, our assumption concerning the Fermi surface 
is considered to be valid for Cul234. 

4. Summary 

We demonstrated that the ground state of the 3D d-p model with a multilayer perovskite 
structure can be in the (i^2_j^2-wave superconducting state up to the second-order in the 
perturbation theory framework. In the multilayer system, the region with large log Age can 
expand further. This is caused by the multilayering effect, which can increase the chance that 
the Fermi surface has VHS and can maintain a high DOS around the Fermi level over a wide 
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n 

Fig. 5. The results for n = 5 and = 0.098. (a) The eigensolution a\{kz). (b) The Fermi surface 
projected onto the plane with kz ~ 0. 

doping region. This multilayering effect works very well when the unit cell contains more 
layers, as long as a well-defined Fermi surface exists. 
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